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Abstract 

In wireless cellular networks with densely deployed base stations, knowing the ve¬ 
locities of mobile devices is a key to avoid call drops and improve the quality of service 
to the user equipments (UEs). A simple and efficient way to estimate a UE’s velocity 
is by counting the number of handovers made by the UE during a predefined time 
window. Indeed, handover-count based mobility state detection has been standardized 
since Long Term Evolution (LTE) Release-8 specifications. The increasing density of 
small cells in wireless networks can help in accurate estimation of velocity and mobility 
state of a UE. In this paper, we model densely deployed small cells using stochastic 
geometry, and then analyze the statistics of the number of handovers as a function of 
UE velocity, small-cell density, and handover count measurement time window. Using 
these statistics, we derive approximations to the Cramer-Rao lower bound (CRLB) 
for the velocity estimate of a UE. Also, we determine a minimum variance unbiased 
(MVU) velocity estimator whose variance tightly matches with the CRLB. Using this 
velocity estimator, we formulate the problem of detecting the mobility state of a UE 
as low, medium, or high-mobility, as in LTE specifications. Subsequently, we derive 
the probability of correctly detecting the mobility state of a UE. Finally, we evaluate 
the accuracy of the velocity estimator under more realistic scenarios such as clustered 
deployment of small cells, random way point (RWP) mobility model for UEs, and 
variable UE velocity. Our analysis shows that the accuracy of velocity estimation and 
mobility state detection increases with increasing small cell density and with increasing 
handover count measurement time window. 

keywords: Cramer-Rao lower bound (CRLB), heterogeneous networks (Het¬ 
Nets), long term evolution (LTE), mobility state estimation, mobile velocity 
estimation, phantom cell, small cells. 
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1 Introduction 


Since the introduction of advanced mobile devices with data-intensive applications, cellular 
networks are witnessing rapidly increasing data traffic demands from mobile users. To keep 
up with the increasing traffic demands, cellular networks are being transformed into het¬ 
erogeneous networks (HetNets) by the deployment of small cells (picocells, femtocells, etc) 
over the existing macrocells. Cisco has recently predicted an 11-fold increase in the global 
mobile data traffic between 2013 and 2018 [^, while Qualcomm has forecasted an astounding 
lOOOx increase in mobile data traffic in the near future [^. Addressing this challenge will 
lead to extreme densihcation of the small cells which will give rise to hyper-dense HetNets 
(HDHNs). 


Mobility management in cellular networks is an important task that is critical to provide 
good quality of service to the mobile users by minimizing the handover failures. Velocity of a 
user equipment (UE) plays critical role on the handover performance of the UE particularly 
when the cell density is high [3]-[^, where knowing the UE’s velocity becomes necessary 
for effective mobility management. In homogeneous networks that only have macrocells, 
handovers are typically hnalized at the cell edge due to large cell sizes. With the deployment 
of small cell base stations (SBSs) [6]-[^, due to smaller cell sizes, it becomes difficult to 
hnalize the handover process at the cell edge for mobile devices 


10 . In particular. 


high-mobility devices may run deep inside the coverage areas of small cells before hnalizing 
a handover, thus incurring handover failure due to degraded signal to interference plus noise 
ratio (SINK). These challenges motivate the need for UE-speci£c and cell-specihc handover 
parameter optimization, which typically require estimation of the UE’s velocity for effective 


conhguration of handover parameters 11 . A UE’s velocity estimate may also be used for 


scheduling 12,13, mobility load balancing 14 , channel quality indicator (CQI) feedback 


enhancements 15 , and energy efficiency enhancements 16 


In this paper, we introduce a novel and efficient handover-coun^ based UE velocity 


^Subsequently, “handover-count” is used in a broad sense to refer to the number of cells traversed by a 
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estimation technique using the tools from stochastic geometry, and characterize its accuracy 
through Cramer-Rao Lower Bounds (CRLBs), when the density of SBSs is known. Since the 
service provider has the information of the number of SBSs in a particular geographic area, 
the SBS density in that area can be calculated and broadcasted as part of system information 
in next generation networks. The SBS density may also be signaled in a user-specihc manner 
to the next generation UEs which are capable of velocity estimation. Our contributions in this 
paper are as follows: 1) for a given small-cell density, two approximations to the probability 
mass function (PMF) of handover-count of a UE are derived using a heuristic approach; 2) 
using the PMF approximations, expressions for the CRLB of velocity estimation are derived; 
3) a minimum variance unbiased (MVU) velocity estimator is derived whose variance tightly 
matches with the CRLB, and accuracy of the estimator is investigated for various UE speeds, 
SBS densities, and handover-count measurement times; 4) the estimated velocity is used to 
detect the mobility state (low/medium/high) of a UE, and the expressions for the probability 
of detection and probability of false alarm are derived; 5) accuracy of the velocity estimator is 
analyzed for different realistic scenarios such as: RWP mobility model for the UEs, clustered 
deployment of SBSs, and variable UE velocity. 

This paper is organized as follows. We briefly review the related works and the state of 
the art in Section]^ In Sectionwe describe our system model of small-cells using stochastic 
geometry. Subsequently, we derive the approximations to the PMF of handover-count of a 
UE in Section In Section we hnd CRLB for the velocity estimation of UE and also 
derive a MVU velocity estimator. In Section]^ we provide expressions for the mobility state 
(low, medium, high) probabilities, the probability of detection, and the probability of false 
alarm. We show numerical results on our Endings in Section and conclude the paper in 
Section |8l 


UE. E.g., in phantom cells 17 , a UE is always connected to the MBSs, and can get additional throughput 


from SBSs, minimizing handover failures. 
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2 Review of Prior Work 


Existing long term evolution (LTE) and LTE-Advanced technologies are capable of es¬ 
timating the mobility state of a UE into three broad classes: low, medium, and high- 
mobility [^|^[^[^. This is achieved at a UE by counting the number of handovers within 
a given time window, and comparing it with a threshold which can be implemented during 


the connected mode 19 or the idle mode 20 of the UE. It can also be implemented at the 
network side, by tracking the prior history of handovers for a particular UE. The accuracy 
of a UE’s mobility state estimate will benefit signihcantly from the densihcation of SBS 
deployment. While more accurate UE-side speed estimation techniques based on Doppler 


estimation have been discussed in 21-24 , due to their complexity and standardization chal¬ 
lenges, they have not been adopted in existing cellular network standards. While global 
positioning system (GPS) can be used for accurate estimation of a UE’s speed, it may not 
be a practical solution for mobility management as 1) the GPS receiver at a UE consumes 
signihcant amount of power, 2) GPS coverage may not be available in environments such as 
urban canyons and underground subways, and 3) not all the UEs are equipped with GPS 
receivers. A popular usecase example that uses such real-time mobility state information at 
a UE is that of UE-speci£c cell selection: high speed UEs can be biased to stay connected 


to macrocells even when small cell link quality is better 25 


There are various studies available in the literature on handover-count based mobility 


state detection (MSD). A simulation based mobility analysis is performed in 18 where it 
is proposed that the number of handovers made by a UE are weighted differently for macro 
to macro, macro to pico, pico to macro, and pico to pico handovers (1, 0.45, 0.25, and 0.1, 
respectively), to produce a good estimate of UE mobility. In j^, considering a random 
way point (RWP) model and using stochastic geometry, the expected number of handovers 
during the movement period of a RWP model is derived. The probability density function 


(PDF) of the sojourn time is also presented for homogeneous networks. In 27 , a theoretical 
framework using stochastic geometry is developed to study the UE mobility in HetNets, in 
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which the expressions for vertical and horizontal handoff rates are derived. In 28 , closed 


form expressions of cross-tier handover rate and sojonrn time in a small cell are provided 


nsing stochastic geometry for a two-tier network. In another stndy 29 , stochastic geometry 
is nsed to derive the handover rate, which in tnrn is nsed to derive coverage probability of a 
UE in HetNet by considering that the UE is mobile and a fraction of the handovers resnlt 
into failnre. 

A set of related works in the literatnre have stndied the path prediction of UEs and the 
handoff time estimation along the predicted path which can help in improving the qnality 
of service to the UEs. A destination and mobility path prediction model called DAMP is 
proposed in [30|, while a handoff time window estimation method and mobility-prediction- 


aware bandwidth (MPBR) allocation scheme is presented in 31 . In 32 , different methods 
for predicting the fntnre location of a UE based on prior knowledge of the UE’s mobility 
are stndied. A method for tracing UE’s location nsing semi-snpervised graph Laplacian 


approach is proposed in 33 


Despite the earlier work in 3 4,11,18-24,26,27,29-33 , fnndamental performance bonnds 


and optimnm algorithms for handover-connt based velocity estimation have not been stndied 
in the literatnre. The focns of this paper is on velocity estimation and MSD based on 
the handover connts of a UE; path prediction and handoff time estimation techniqnes are 
not considered. Mobility state estimation itself is an important topic in 3GPP Release-8 
specifications, and being researched cnrrently. Additionally, there are applications which 
reqnire velocity estimation at the UE side, in which case, implementing the path prediction 
algorithms at a UE would be difficult. To our best knowledge, velocity estimation based on 
handover counts at a UE has not been studied analytically in the literature, which is one of 
the main contribution of this paper. 
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3 System Model 


Consider a HetNet scenario where the macro base stations (MBSs) tier and the SBSs tier nse 
different freqnency bands, as in phantom cell architectnre (^. We assnme that the SBSs are 
randomly distribnted according to a homogeneous Poisson point process (PPP) with intensity 
A. As the SBSs use a dedicated frequency band which is different from the frequency band 
used by the MBSs, the coverage of a small cell would depend on the neighboring SBSs, 
but not on MBSs. Hence, the small cells of the network can be modeled using Poisson- 
Voronoi tessellation as illustrated in Fig. [^b). In Fig. [^a), co-channel deployment of small 
cells is illustrated in which the coverage of small cells is determined by the interference 
from neighboring SBSs and MBSs. In this paper, we consider only the dedicated channel 
deployment of small cells as illustrated in Fig. [^b). 




Figure 1: Illustration of cellular network layout; (a) coverage of macrocells and small cells 
in co-channel deployment, as in a typical HetNet; (b) coverage of small-cells with dedicated 
channel deployment as in a phantom cell network. 

We consider a simple UE mobility scenario as also shown in Fig. Bb) , in which the 

UE travels along a linear trajectory (for example, through X-axis) with constant velocity v. 

During the travel, we assume that the UE can determine whenever it crosses the boundary of 

a small cell. In a broad sense, we call the boundary crossings made by the UE as handovers. 

Therefore, the number of handovers H made by the UE during a measurement time window 

T is equal to the number of intersections between the UE travel trajectory (of length d = vT) 
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and the small-cell boundaries. We use linear mobility model for its simplicity in theoretical 
analysis, and this model is suitable for scenarios such as medium/high speed cars and trains 
that may travel through downtown areas. There may be many small cells deployed in such 


urban areas in the future which may be referred to as ultra/hyper-dense networks 34 . Linear 
mobility may not be accurate for some other scenarios, therefore we have also considered 
RWP mobility model which is more general and includes linear mobility as a special case. 

Note that the approach in Fig. [^b) for handover-count based MSD of a UE, into 
low/medium/high mobility states, has already been specihed in the LTE Release-8 stan¬ 
dard. Remarkably, no studies exist in the literature that investigate fundamental bounds 
and effective estimators for UE velocity. At high UE speeds and high SBS densities, han¬ 
dover failure of a UE becomes more likely [5,35,36 . Hence, the velocity estimation of UE 
based on handover counts may not work effectively with conventional LTE mobility man¬ 


agement. On the other hand, emerging small-cell architectures such as phantom cells 17,37 


decouple the control and user planes, and allow the UE to be connected to macrocell all the 


time. Small-cells are discovered through special discovery signals 38 , and a UE can connect 
them (reminiscent to a handover) to have higher throughput. Hence, handover-count based 
velocity estimation is still applicable in such scenarios. 

The handover count H for a scenario as in Fig. [^b) is a discrete random variable and 
its statistics do not change with the direction of the linear trajectory because the SBS loca¬ 
tions are modeled using a homogeneous and stationary PPP. The probability mass function 
(PMF) fuih) of handover count obtained through the simulations is shown in Fig. |^a)-|^d) 
for different SBS density A and UE velocity v settings, with the handover-count measure¬ 
ment time hxed to T = 12 s. For low SBS densities, as in Fig. [^a), the PMFs for different 
UE velocities are overlapping signihcantly, leading to low velocity estimation accuracies. For 
higher SBS densities, as in Fig. j^d), the PMFs for different velocities are separated, lead¬ 
ing to better estimation accuracies. It should also be noticed that for higher UE velocities, 
standard deviation of the handover count also increases, implying lower estimation accura- 
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(a) 



(c) 



(b) 



(d) 


Figure 2: PMF of the handover count for different A and v values; (a) A = 100 SBSs/km^; 
(b) A = 200 SBSs/km2; (c) A = 500 SBSs/km^; (d) A = 1000 SBSs/km^. 


cies. Another characteristic of the PMF is that for large A and v values, the shape of PMF 
resembles Gaussian distribution. 


3.1 Modeling Handover-Count Statistics 


In order to obtain the velocity estimate of a UE based on its handover count H, we need to 
know the PMF fnih) of handover count. For a scenario as in Fig. [^b), an exact expression 


for the mean number of handovers can be derived as 26,39 




4vT^ 


TT 


( 1 ) 


To the best of our knowledge, there is no expression available for the PMF of handover count 
in the literature. Deriving an expression for the PMF fnih) is a complicated and laborious 


















task, which might result into a mathematically intractable expression 39 . Hence, we derive 
an approximation to the PMF fnih) in this paper. 

There are several papers in the literature where the PPP based parameters are approx¬ 
imated rather than deriving the exact expressions due to the complexity involved in the 


derivation of exact expressions. For example, in 40 , geometrical characteristics of the 
perimeter, area and number of edges in a 2-dimensional Voronoi cell, and volume, surface 
area and number of faces in a 3-dimensional Voronoi cell are approximated by htting gen¬ 
eralized gamma distribution to the respective histograms. Similarly, the distributions of 


2-dimensional cell area and 3-dimensional cell volume are approximated in 41 by htting 
some simple expressions. 

In this paper, we derive two approximations to the PMF fnih) for the handover count: 

1) approximation ffj{h) derived using gamma distribution; 

2) approximation f^ih) derived using Gaussian distribution. 

These two approximations to the handover count PMF will be discussed in more detail in 
Section]^ and their accuracies will be further investigated and compared in Section 7.1 


4 Approximation of the Handover Count PMF Using 
Gamma and Gaussian Distributions 

In this section, two approximations for the PMF of handover count will be introduced. 
In each approximation method, the parameters of a distribution (gamma distribution or 
Gaussian distribution) will be approximated using the curve htting tools in Matlab. 
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4.1 Approximation of the PMF of Handover Count using Gamma 


Distribution 


Gamma distribution has been commonly used in approximating the statistical distribution 
of the parameters related to PPPs such as area, volume, number of edges, etc., of Poisson 


Voronoi cells 40,41 . It can be effectively used to approximate the handover count PMF. 
The gamma PDF can be expressed using the shape parameter a > 0 and rate parameter 
/? > 0 as 


p{x) = for x € (0,cx)), (2) 

F(a) 

where F(a) = is the gamma function. Fitting the gamma PDF to the PMF 

fnih) is not a straight forward task since the gamma PDF is a continuous function while 
the PMF fnih) is a discrete function. Therefore, we £t the subsets of gamma PDF to the 
PMF fnih). The subsets of gamma PDF can be obtained by integrating the gamma PDF 



Figure 3: Illustration of htting gamma distribution to the handover count PMF. 
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between the integer values of x as illustrated in Fig. Integrating the gamma PDF between 
0 and 1 provides the PMF value for h = 0, integrating the gamma PDF between 1 and 2 
provides the PMF value for h = 1, and so on. This process can be mathematically described 
as: 



f^{x)dx, for h e {0,1, 2,...}, 


( 3 ) 


where, f%{h) is an approximation to the PMF fnih). Substituting (|^ into ([^, we get 


fUh) = 


Qa 




T{a 


-X 




/S' 


a nh-\-l 


F(«) 


for h e {0,1, 2,...}. (4) 


With the change of variable x = t/jS, we can rewrite Q in its equivalent form as 


fUh) 


pa i>l3{h+l) 

r(a) 7/3/1 



r(«) V F(a) 


( 5 ) 


where, F (a, (3h, (3{h + 1)) = is the generalized incomplete gamma function. 

However, for the approximation in (|^ to be accurate, the values for a and j3 parameters 
should be chosen such that the mean squared error (MSE) between fnih) and f%{h) is 
minimized. 


Lemma 1. The a and [3 parameters for the approximation in (|^ which minimize the MSE 
between fnih) and ffj{h) can he expressed as 


2.7 + 4dyA, (6) 

0-8 

TT H- -J=^ (7) 

0.38 Td^A 

where d = vT is the distance traveled by UE during the handover-count measurement time. 
Proof. See Appendix □ 
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Using (§- 0 . we can capture the statistical distribution of handover count if the SBS 
density and the distance traveled by UE are known. Through ([^ and Q, it can also be 
noticed that the handover-count distribution depends on the distance traveled by UE, rather 
than the UE velocity or the time period independently. 

4.2 Approximation of the PMF of Handover Count using Gaus¬ 
sian Distribution 

In the previous section, approximation to the handover count PMF fnih) was derived using 
gamma distribution which resulted into an expression in integral form. In this section, we 
will approximate the PMF fnih) using Gaussian distribution which results into a closed 
form expression for the handover-count PMF. 

The PDF of Gaussian distribution can be expressed as a function of mean /x and variance 


r(^) 


1 _ 

P. 2 < t ^ 

\/27ra‘^ 


( 8 ) 


Since the PMF fnih) is discrete and the Gaussian distribution in (|^ is continuous, we 
consider only non-negative integer samples of the Gaussian distribution for the htting process. 
Henceforth, the approximation to fnih) can be expressed as. 


fnih) = for h e {0,1, 2,...}. (9) 

v27rcr"' 

The values of /i and should be chosen to minimize the MSE between fnih) and /]f(h). 
Lemma 2. The jj, and parameters for the approximation in which minimize the MSE 
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between fnih) and ff^ih) can he expressed as 

Ad^/X 

li , 

TT 

= 0.07 + 0.41ci\/A, 

where d = vT is the distance traveled by the UE during the handover count measurement 
time. 

Proof. See Appendix □ 


( 10 ) 

( 11 ) 


5 Cramer-Rao Lower Bound for Velocity Estimation 


CRLB can be used to serve as a lower bound on the variance of an unbiased estimator 42 


An estimator is said to be unbiased if its expected value is same as the true value of the 
parameter being estimated. An unbiased estimator whose variance can achieve the CRLB is 
said to be an efficient estimator, and it can achieve minimum MSE among all the unbiased 
estimators. In some scenarios, it might not be feasible to determine an efficient estimator. In 
that case, the unbiased estimator with the smallest variance is said to be a MVU estimator. 

In this section, we will derive the CRLB for velocity estimation using the mathematical 
tools from estimation theory. Since we have two different approximations for the PMF of 
the number of handovers, we will obtain two separate CRLB expressions. 


5.1 CRLB Derivation using Gamma PMF Approximation 


In this sub-section we will obtain the CRLB by considering the handover count PMF ap¬ 


proximation that was derived using gamma distribution in Section 4.1 


Theorem 1. In a Poisson-Voronoi tessellation of small cells with SBS node density X, let a 
UE travel with velocity v over a linear trajectory and make H handovers over a time duration 
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T. If the PMF of the handover count can he expressed using ffj{H;v) as in ([^, then the 
CRLB for velocity estimation is given by 


var('y) > 


E 


/ aiogffj(H-,v) 
dv 


( 12 ) 


where, E[-] is the expectation operator with respect to H, and 


d\ogfUf^;v) 

dv 




ck; ck + 1, a + 1; —/dh) 


0(2 r(Q(, (3h,(3{h + 1)) 

— {h -\- l)'^2-f2(Q(, o; o + 1, o + 1; — ld(h + 1)) 
AT^^\ 


r(Q(, (3h, /3{h + 1)) 


'y{a,/3h) \og{/3h) - -f{a, /3{h + l))\og{/3{h + 1)) 


0.8TVX/3°-^e-f^^ - e-f^{h + 1)“1 ^ ^ ^ 

-^-4TvA fjia), 

T{a,f3h,f3{h + l))i0.38 + vTV\y 


(13) 


where, is digamma function,{a, x) = Jq t^" ^dt is lower incomplete gamma function, 
and 2 F 2 {ai,a 2 ]hi,h 2 ]z) is generalized hypergeometric function which is expressed as 


O' 2 'i hi,h 2 , z) — 


^ {bi)k{b2)k k\' 


(14) 


where, (a)o = 1 and {a)k = a{a + l)(a + 2)...(a + k — 1), for k > 1. 


Proof. See Appendix 


□ 


Due to the complexity of expression in (13), it is impractical to derive the right hand 


side (RHS) of (12) in closed form. For this reason, we can only hnd asymptotic CRLB by 


numerically evaluating the RHS of (12). Through simulations, we generate N samples of the 


random variable H and denote them as {Hn}, for n G 1,2,...,A^. Using these N samples. 
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we can numerically evaluate the asymptotic CRLB using 


var(t)) > 


N 


( Af f 9logfH(m-,v) 
Z^m = Hmin I V Sv 



( 15 ) 


where, hTmax = max{if„ : Vn G 1, 2,A^}, is the maximum value of -ffmin = niin{iif„ : 

Vn G 1, 2,A^}, is the minimum value of Hn, and N^, = '^n=i is the number 

of elements in the set {Hn} that are equal to m. Here, !{■} is the indicator function whose 
value is 1 if the condition inside the braces is true, 0 otherwise. 


5.2 CRLB Derivation using Gaussian PMF Approximation 

In this sub-section, we will obtain the CRLB by considering the PMF approximation using 
Gaussian distribution that was derived in Section 14.21 

Theorem 2. In a Poisson-Voronoi tessellation of small cells with SBS node density X, let a 
UE travel with velocity v over a linear trajectory and make H handovers over a time duration 
T. If the PMF of the handover count can he expressed using f}j{H;v) as in 0 , then the 
CRLB for velocity estimation is given by 


var('i)) > 




(16) 


Proof. Consider the PMF approximation ffjih) in ([^ which can be represented as a general 
Gaussian distribution, H J\f (jj, a^) , where /i and are given by (10) and ( [II| respectively. 
The Fisher information for the general Gaussian observations is given by Section 3.9] 




dv J 2 (cr^)^ y 




iTy/X \ 1 1 


TT 


cr" 2 (a^ 


(^0.41TyA)^ (17) 


Using inverse of the Fisher information, the CRLB can be expressed as in (16). 


□ 
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5.3 Minimum Variance Unbiased Estimator for UE Velocity 


In Section 4.1 and Section 4.2 two CRLB expressions were derived by considering gamma 
and Gaussian distributions, respectively, for approximating the handover count PMF. In the 
case of using gamma distribution, the CRLB expression was complicated and not in closed 
form. On the other hand, in the case of using Gaussian distribution, the CRLB expression 
was relatively simple and in closed form. Hence, in this sub-section, we will consider the 
case with Gaussian distribution and derive an estimator v for a UE’s velocity, which takes 
the number of handovers H as the input. We will further derive the mean and the variance 
of this estimator and show that it is a MVU estimator. 

To derive the MVU velocity estimator, we hrst use Neyman-Fisher factorization to hnd 


the sufficient statistic for v 42, Section 5.4]. Then, we make use of Rao-Blackwell-Lehmann- 
Scheffe (RBLS) theorem to hnd the MVUE Section 5.5]. The Neyman-Fisher factoriza¬ 
tion theorem states that if we can factor the PMF f^ih) as 


JhW = g(T{h),v)r(h), 


(18) 


where g is a. function depending on h only through J^{h) and r is a function depending only 
on h, then J^{h) is a sufficient statistic for v. Using (|^ and letting J^{h) = h, we can factor 


the PMF ffiih) in the form of (18) as 


fnih) = ^=e 2 -" 1 


72 




(19) 


g{T(h),v) r(h) 


Therefore, the sufficient statistic for v is J^{h) = h. The sufficient statistic can be used 
to hnd the MVU estimator by determining a function s so that v = s{J^) is an unbiased 
estimator of v. By inspecting the relationship between the mean number of handovers H 
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and the velocity n in Q, we can formnlate an estimator for v as: 


V = 


ttH 

AT^/X' 


(20) 


In order to evalnate whether this estimator is nnbiased, the expectation of the above esti¬ 
mator can be derived as 


E[v] = E 


ttH 


AT^/X\ ATy/X 


"" m] = ^E 


at^A 


( 21 ) 


Plngging (10) into (21), we get 


P TT AvTA 

E\v\ = -=-= V. 

AT A tt 


( 22 ) 


Therefore, the estimator v expressed in (20) is nnbiased. Since this estimator is derived 


throngh RBLS theorem, it is an MVU estimator. To determine whether it is an efficient 
estimator, we derive the variance of the MVU estimator as follows: 


var(n) =var -= 

\atAJ 




(23) 


Comparing (23) with (16), we can notice that the variance of MVU estimator is greater than 


the CRLB, and hence, the derived estimator is not an efficient estimator. Nevertheless, in 


Section 7.3, we show that the variance of the MVU estimator is very close to the CRLB. 


6 Mobility State Detection 


In this section, we will perform statistical analysis of MSD, in which a UE is categorized into 
one of the three different mobility states: low mobility, medinm mobility and high mobility. 


as in 3GPP LTE Release-8 specihcations 18,43,44 . We assnme that the nnbiased velocity 


estimator derived in Section 5.3 is used to estimate the UE velocity, and we will derive 
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expressions for the probabilities that a UE is categorized into each of the three mobility 
states. 


Using the estimated UE velocity i) from (20), the UE can be categorized into one of 
the three mobility states: low (Sl), medium (S'm), and high (Su), based on the following 
conditions: 


5 


5'l if h < vi, 

< if vi < i) < Uu, 

Sn if h > v^, 


(24) 


where, S G {S'l, Sm, 5'h} is the detected mobility state of the UE. The thresholds ui and 
are the lower and upper velocity thresholds, respectively, based on which a UE is classihed 
into one of the three mobility states. 


6.1 Mobility State Probabilities 

For a given velocity v, we dehne mobility state probability as the probability that the UE 
is categorized into a particular state. We can dehne the following three mobility state 
probabilities: 


P{S = Si^;v) —^-Probability that the mobility state is detected as S'l, for a velocity v; 

P{S = S'm; v) —)-Probability that the mobility state is detected as Sm, for a velocity v; 

P{S = S'h; v) —)-Probability that the mobility state is detected as Sh, for a velocity v. 


For a given velocity v, as the number of handovers H is a random variable, the velocity v 


estimated using (20) is also a random variable. Hence, there will be false alarms and missed 
detections for calculating the mobility state. 


Next, we derive analytic expressions for the mobility state probabilities. Using (20), we 
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can express the PMF of v as, 


h{v) =P{v = v)=P 


f irH \ ^ ATy/Xu \ , f AT^Xv 


\AT^^X 


= iy \ =P \ H = 


TT 


— fn 


v ^ 


= fH{h), (25) 


where, h — Using the approximation of PMF with Gaussian distribution as i 

(§. we can approximate the PMF of v as, 


fv{l^) =fH{h) 

^fnih) 


V2 


1 

e 2 <t^ , for h G {0,1, 2,...}. 




(26) 


Now, we can express the three mobility state probabilities as 

hi 


P(S = = P(v < V,) K 5^/S(A), 

h=0 

P(S = Sm;v) = P(vi<v<Vi,)^ f^(h), 

h=hi-i-l 

oo 

P(S = S„;t>) = P(f.>t>„)« 5^ /S(h), 

where, hi = 


h=hu-hl 


(27) 

(28) 

(29) 

(30) 


are the optimnm lower and upper handover count thresholds for MSD, respectively. Given 
the velocity thresholds vi and Uu, the choice of handover count thresholds has a direct im¬ 
pact on the probability of correctly detecting the mobility state of a UE based on its ve¬ 


4T^ni 


AT^Avh 

TT 

j CXilGl / 

TT 


locity 45 . In (30), we have theoretically derived the handover-count thresholds for MSD 


which are optimum for the given velocity thresholds. In other related works in the liter¬ 
ature, the handover-count thresholds for MSD have been determined through simulations. 
In [5,45,46 , the handover-count thresholds for MSD are found heuristically by considering 
the cumulative distribution function (CDF) plots of handover counts, for few different UE 


velocities. In 47 , the handover counts are assumed to be distributed as Gaussian PDF, and 
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the optimum handover-count thresholds are obtained by using the handover-count PDFs 
for few different UE velocities. However, in these prior works, the optimum handover-count 
thresholds are determined for some particular values of BS density and measurement time 
windows. Moreover, the statistical relationship between the UE velocity and the handover 
count is not considered. In this paper, we have derived general expressions for the optimum 
handover-count thresholds as a function of SBS density A, handover-count measurement time 
T, and velocity thresholds (ui and Uu). 


6.2 Probability of Detection and Probability of False Alarm 


The probability of detection is the probability that the mobility state of a UE is detected 
correctly. Mathematically, it can be expressed as 


Pn = 


P{S = Si^;v) if V < v\, 

P{S = S'm; v) if v\<v< Uu, 


P{S = Syi]v) if n > Uu- 


(31) 


The probability of false alarm is the probability that the mobility state is detected incorrectly, 
which can be mathematically expressed as PpA = 1 — Pd- 

Consider an illustrative example in which the UE velocity is n = 60 km/h, SBS density 
is A = 1000 SBSs/km^, and handover-count measurement time is T = 12 s. The PMFs 
fnih) and fv{h') are shown in Figs. I^a) and|^b), respectively, which were obtained through 
Monte Carlo Simulations. In Fig. |^b), the range of v is divided into three regions S'l, Sm 
and Fh that are separated from each other through the velocity thresholds v\ = 40 km/h 
and Uu = 80 km/h. It can be noticed in Fig. |^b) that even though the actual UE velocity 
is a constant u = 60 km/h which belongs to S'm state, the estimated velocity is spread into 
a range of velocities. Hence, there is a small probability that the mobility state could be 
erroneously detected as S'l or S'h states, which is the probability of false alarm PpA. On 
the other hand, majority of the times the mobility state would be correctly detected as Sm, 
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which is the probability of detection Pd- 




No. of handovers, H Estimated UE velocity, v 

(a) (b) 

Figure 4: (a) PMF of the number of handovers; (b) PMF of the estimated velocity; v = 
60 km/h, A = 1000 SBSs/km^, T = 12 s, ui = 40 km/h and Vu = 80 km/h. 


With the particular parameter settings shown in Fig. the probabilities of different 
mobility states can be evaluated using (27)-(|2^ as P{S = 5l;u = 60) = 0.047, P{S = 


Sm',v = 60) = 0.8821, P{S = S'h;u = 60) = 0.0709, respectively. Then, the probabilities 
of detection and false alarm can be evaluated using ( [3l| as Pd = P{S = Sm',v = 60) = 
0.8821, Pfa = P(5 = 5'l; v = 60) + P{S = S'h; v = 60) = 0.1179, respectively. 


7 Numerical Results 


In this section, firstly, we will validate the accuracy of gamma PMF approximation ffj{h) 
and Gaussian PMF approximation ffi{h) by plotting their MSE performances. Secondly, 
we will plot the CRLBs and analyze the achievable accuracy of a UE’s velocity estimate for 
different SBS density A, UE velocity v, and handover count measurement time T. Finally, 


we will plot the variance of the MVU velocity estimator derived in Section 5.3 and show that 
it is approximately equal to the CRLB. 


7.1 Accuracy of PMF Approximation 

Approximations to the handover count PMF were derived using gamma distribution and 


Gaussian distribution in Sections 4.1 and 4.2, respectively. In this sub-section, we will 
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Figure 5: MSE versus A for different u, and T = 12 s; (a) with PMF approximation using 
gamma distribution; (b) with PMF approximation using Gaussian distribution. 

quantify the accuracy of each approximation method by evaluating the MSE between the 
approximation and the PMF fnih). The MSE can be expressed as 



(32) 


where iVh is the number of points in the PMF. The characteristics of MSE with respect to 


the variations in A and v are shown in Fig. for the two approximation methods. In gen¬ 


eral, the MSEs of both approximation methods decrease with the increase in SBS density 
A or UE velocity v. In other words, higher SBS density and higher UE velocity leads to 
better accuracy of the PMF approximation. By comparing Fig. |^a) and Fig. |^b), it can be 
noticed that the approximation using gamma distribution provides approximately ten times 
smaller MSE than the approximation using Gaussian distribution. However, the approxima¬ 
tion using gamma distribution is not in closed form and hence it is more complicated than 
the approximation using Gaussian distribution. Therefore, there exists a trade-off between 
the accuracy and the complexity while making a choice between the two approximation 
approaches. 
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Figure 6: CRLB versus v for different A, and T = 12 s; (a) CRLB derived using f%{h)] (b) 
CRLB derived using f^ih). 


7.2 CRLB Results 


The CRLB plots for UE velocity estimate v are shown in Fig. for the variations in UE 
velocity v and SBS density A. In the case of gamma approximation of PMF fnih), the 
CRLB plots in Fig. a) are obtained by numerically evaluating the expression in ( [l5| ). On 
the other hand, the CRLB plots in Fig. |^b) are obtained using the closed form expression 


in (16), for the case of Gaussian approximation of PMF /^(h). It can be observed that 


the CRLB plots in Fig. |^a) and Fig. |^b) are similar and follow same trends with respect 
to A and v. However, we can expect that the plots in Fig. [^a) are more accurate because 
of the smaller MSE of gamma approximation method. On the other hand, the closed form 


CRLB expression in (16) for the case of Gaussian approximation can provide more insights. 


for example, on the dependence of GRLB into different parameters. 

From Fig. it can be noticed that the GRLB increases with the increasing UE velocity 
V. This is because the variance of the number of handovers H increases with increasing 
V, which can be observed in Fig. In contrast, the GRLB decreases with increasing SBS 
density A, which can also be intuitively understood from Fig. With A = 100 SBSs/km^, 
the peaks of the PMFs for different UE velocities are closely spaced with each other making it 
difficult to distinguish between the different UE velocities, resulting into higher GRLB. With 
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Figure 7: CRLB versus T for different A and u = 60 km/h; (a) CRLB derived using PMF 
approximation (b) CRLB derived using PMF approximation f^ih). 


A = 1000 SBSs/km^, the peaks of the PMFs have more spacing between them making it 
easier to distinguish between the different UE velocities, resulting into lower CRLB. Results 
in Fig. I^show that for all considered SBS densities, using a handover count based approach, 
a UE’s velocity can be estimated with a root mean squared error (RMSE) of 20 km/h for 
UE velocities less than 40 km/h. For UE velocities on the order of 120 km/h, velocity can 
still be estimated with a RMSE less than 31 km/h for densities of 100 SBSs/km^, which can 
be further improved for higher SBS densities. 


In related works of the literature 18,45-47 , a hxed measurement time interval of 30 s, 
60 s, or 120 s is used for the MSD of UEs. However, we have used a smaller measurement 
time interval of T = 12 s in Fig. [^so that the velocity estimator can provide quicker results. 
At the same time, T = 12 s also provides a reasonable estimation accuracy. For example, 
with SBS density A = 500 SBSs/km^ and a UE velocity v = 120 km/h, the CRLB of velocity 
estimation is just 20 km/h as shown in Fig. After all, it is a choice of the service provider 
to use a measurement time interval based on the requirements. Therefore, the effects of 
variations in T are also investigated in the remainder of this paper. 

The effect of handover-count measurement time T on the CRLB can be seen in Fig 
For a given UE velocity and SBS density, CRLB decreases as T increases. Therefore, longer 
handover-count measurement time results into better accuracy of velocity estimation, with 
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the assumption that the UE will continue traveling on a linear trajectory. 


7.3 Variance of the MVU Velocity Estimator 


The variance of the MVU velocity estimator given in (23) is plotted in Fig. Even though 
the derived MVU estimator is not an efficient estimator, it can be observed from the plots 
that the variance of the MVU estimator tightly matches with the CRLB. Therefore, the 
MVU estimator performs very close to an efficient estimator. 




Handover count measurement time, T (s) 

(b) 

Figure 8: Comparison between the CRLB and the variance of the MVU estimator; (a) for 
variations in v and A, with T = 12 s; (b) for variations in T and A, with u = 60 km/h. 


7.4 Mobility State Probabilities, and Probabilities of Detection 
and False Alarm 

In order to classify a UE into one of the three mobility states, the velocity thresholds v\ and 
Uu can be set depending on the requirements of the service provider. In this paper, we have 
set Ui = 40 km/h and v-^ = 80 km/h. The plots of mobility state probabilities versus UE 
velocity u, for four different A values are shown in Fig. Here, the analytic plots that were 


obtained using (27)-(29) are observed to match accurately with the simulation results. When 


we examine these plots carefully, it can be noticed that as the SBS density A increases, the 
slope of the curves during their transitions also increases, which is closer to the ideal case 
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Figure 9: Mobility state probabilities versus UE velocity; T = 12 s, ui = 40 km/li and 
Uu = 80 km/h; (a) A = 100 SBSs/km^; (b) A = 200 SBSs/km^; (c) A = 500 SBSs/km^; (d) 
A = 1000 SBSs/km2. 


where the slope is equal to inhnity. In other words, accuracy of MSD improves signihcantly 
with larger SBS densities. 


The probabilities of detection and false alarm versus the UE velocity v are shown in Fig. 10 


for different A values. As the UE’s velocity nears to any one of the velocity thresholds, 
the probability of detection decreases and the probability of false alarm increases. The 
probability of detection is close to 0.5 when the UE velocity is equal to ui or Uh- 


The average probability of detection is plotted in Fig. H for different combinations of the 
handover count thresholds hi and such that h^> hi. Here, the probability of detection is 
averaged over UE velocity in the range of 10 km/h to 120 km/h, with A = 500 SBSs/km^. It 
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Figure 10: (a) Probability of detection; (b) probability of false alarm. 


can be observed in Fig. 11 that the average probability of detection maximizes to 0.797 with 
hi = 3 and = 7. On the other hand, computing (30) with A = 500 SBSs/km^, T = 12 s, 
v\ = 40 km/h and Uu = 80 km/h also results into hi = 3 and hu = 7. Therefore, equations in 


(30) obtained as a result of our analysis provide analytic expressions for optimum h\ and h„ 
that maximizes the probability of detection. 



Figure 11: Average probability of detection for different combinations of hi and h^, with 
A = 500 SBSs/kmT 


7.5 Estimator Performance with RWP Mobility Model 


In this sub-section, we assess the accuracy of velocity estimator with RWP mobility model 


for the UE 26 . In the RWP model, the UE movement trace is assumed to have a sequence 
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where n denotes the n-th 


of quadruples which can be defined as 
movement period, Xn-i and Xn denote the starting and target waypoints, respectively, dur¬ 
ing the n-th movement period, 14 denotes the velocity, and Sn denotes the pause time at 
the waypoint The angle between two consecutive waypoints is uniformly randomly dis¬ 
tributed on [—TT, tt], while the transition length = ||X„ — X„_i|| between two consecutive 
waypoints is i.i.d. and Rayleigh distributed with the CDF, P{L <l) = 1 —exp (—I > 0, 
where ^ is dehned as the mobility parameter. Larger ^ statistically implies that the transition 
lengths L are shorter and may be appropriate for mobile users walking. In contrast, smaller 
^ statistically implies longer transitions lengths which may be appropriate for driving users. 

We performed simulations for a special case of the RWP mobility model in which the UE 
velocity Vn = v is a. positive constant, and the pause times Sn = 0. The characteristics of 


the RMSE of velocity estimator are shown in Fig. 12 It can be observed that the RMSE 
increases with increasing v, and decreases with increasing T, similar to the characteristics 
of the linear mobility model. The RMSE increases with the increasing mobility parameter 
4 For large 4 the UE switches its directions more number of times within the handover 
count measurement time interval, leading to larger estimation errors. On the other hand, 
for smaller 4 the UE’s direction switch rate is smaller which results into smaller RMSE. As 
^ 0, the RMSE of RWP mobility model converges to the RMSE of the linear mobility 

model. This is because the direction switch rate tends to zero, and the UE follows a straight 
line trajectory indehnitely. 


7.6 Mobility State Detection with Variable UE Velocity 

In this sub-section, we demonstrate the functionality of our MVU estimator with variable UE 
velocity, and the effect of handover-count measurement time on estimation accuracy. Con¬ 
sider an example in which a user is traveling in a train that is moving over a straight line tra¬ 
jectory in a downtown region. Assume that the density of small cells is A = 1000 SBSs/km^, 
and handover counts are measured during regular intervals T. The actual velocity of UE 




(a) (b) 

Figure 12: RMSE performance of the velocity estimator with RWP mobility model; (a) 
RMSE versus UE velocity, with T = 12 s and A = 500 SBSs/km^; (b) RMSE versus handover 
count measurement time, with v = 60 km/h and A = 500 SBSs/km^. 


in the train, estimated velocity and the detected mobility state are shown in Fig. 13 during 


a trip from one station to another station. It can be noticed that the accuracy of velocity 


estimation is better in Fig. 13 b) with T = 30 s, when compared to Fig. [I3[a) with T = 12 s. 
Also, the probability of false alarm is smaller with T = 30 s. However, longer measurement 


time intervals can lead to lower estimation accuracies if the velocity changes signihcantly 


within the measurement time interval. 



0 5 10 15 

Time (min.) 
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0 5 10 15 
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Figure 13: Demonstration of velocity estimation and MSD with variable UE velocity; (a) 
T = 12 s; (b) T = 30 s. 
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7.7 Velocity Estimation in Clustered SBS Deployment 


In practical scenarios, SBSs are deployed more in places where more traffic demand is ex¬ 
pected, which results into clustered deployment. In this subsection, we model the clustered 
deployment of SBSs using Matern-cluster process that can be realized using cluster centers 
and cluster points. The cluster centers form a Poisson process with intensity Aq, and the 
cluster points (SBS locations) are located around each cluster center within a disc of radius 
R. The cluster points in these discs are realized using Poisson processes with intensity Ai. 




Figure 14: RAISE of velocity estimation in homogeneous and clustered SBS deployment; 
(a) RAISE versus the velocity for different Aq, with Ai = 500 SBSs/km^, R = 200 m, and 
T = 12 s; (b) RAISE versus the cluster radius for different Ai, with Aq = 50 clusters/km^, 
n = 60 km/h and T = 12 s. 

The RAISES of velocity estimator in homogeneous and clustered deployments of SBSs are 
compared in Fig. For fair comparison between the homogeneous and clustered deploy¬ 
ments, the intensities of both Poisson and Alatern-cluster processes were kept same. Given 
the parameters Aq, Ai and R of Alatern-cluster process, the equivalent homogeneous intensity 
A for Poisson process can be found using A = AqAittR^. In general, the RAISE is higher in 
clustered deployment of SBSs due to the heterogeneous nature of Alatern-cluster process 
and the fact that a UE cannot obtain handover-count measurements outside the discs of 
clusters. It can be observed in Fig. [I^a) that the RAISE difference between homogeneous 
and clustered SBS deployments increases with increasing UE velocity v. The RAISE dif- 
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ference also increases significantly when Aq decreases, for example, from 100 clusters/km^ 


to 10 clusters/km^. In Fig. 14 b), it can be seen that the RMSE difference decreases with 
increasing R. 


7.8 Velocity Estimation with Matern Hardcore Process for SBS 
Locations 

In this sub-section, we model the SBS locations using Matern hardcore process (HCP), in 
which, the distance between any two points is greater than a hardcore distance i?hc- The 
realization of Matern HCP involves generation of points using a homogeneous PPP, followed 
by a thinning procedure. For the Matern HCP of type H, the thinning procedure involves 
associating a random mark to each point of the parent PPP, and a point is deleted if there 


exists another point within the hardcore distance i?hc with a smaller mark 48 . To realize 
the HCP with intensity Ahc, the intensity of the parent PPP should be 


A = 


ln(l - AhcVrR^J 


(33) 


provided that i?hc is smaller than = l/-\/7rAhc- We use i?hc = p! VvrAhc, where 0 < p < 1 
is the randomness parameter. When p = 0, i?hc is 0 and the HCP is equivalent to PPP in 
which the points are random located as shown in Fig. [I5|(a). When, the p increases to 0.99, 
the points are located more regularly as shown in Fig. [I5|(b). 

The RMSE of the MVU velocity estimator against the variation in p is shown in Fig. IR 
for different velocities. The RMSE of the velocity estimator decreases as the p increases, 
because the randomness in the SBS locations decreases with increasing p. Therefore, in 
scenarios where the SBSs are more uniformly placed, the velocity estimation will work more 
effectively. 


31 








• . . • •• 

•• • •• 

• 

• • 


•• ‘ . • . • • . 

* • • • • 

• . *. • • 

• • • • • 


• • • • . • 

> • • • • 

W 

» • , • • 

• • • 

X 

• • . • 

• • • • • 

CO 

• . • • • 


1 

• • • • •. * 

• • . • . • 

• • • • 


. . 

• : . • • • 

• •• • 

• 

• • • • 


• • • • 

- • • • • • 

• •. • . • • 

• • • • 


X-axis X-axis 


(a) (b) 

Figure 15: Realizations of Matern HCP of type II with fixed Ahc and different p values; 
(a) p = 0; (b) p = 0.99. 



Figure 16: RMSE of velocity estimation with type II Matern HCP for the SBS locations. 


8 Conclusion 


In this paper, we derived a handover-count based MVU estimator for the velocity of a UE, and 
derived optimum handover-count thresholds for MSD. We provided approximations to the 
PMF of handover-count using two types of standard distributions: gamma distribution and 
Gaussian distribution. We found that the PMF approximation using gamma distribution, 
even though complicated when compared to Gaussian distribution, provides a better £t 
to the PMF with a smaller MSE. On the other hand, the use of Gaussian distribution 
provides simpler and closed-form solutions to the PMF approximation. Subsequently, we 
derived the CRLB of velocity estimation and a MVU velocity estimator whose variance is 
approximately equal to the CRLB. The CRLB of velocity estimation decreases with the 
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increasing time interval for counting the number of handovers, which shows the trade-off 
between the accuracy and the rapidness of velocity measurements. 

The increasing density of small-cells in the future cellular networks is facilitates more 
accurate UE velocity estimation, since the CRLB decreases with the increasing SBS den¬ 
sity. Moreover, the probability of MSD also increases with the increasing SBS density. The 
results with RWP mobility model showed that the accuracy of the velocity estimator de¬ 
creases with the increasing randomness in UE’s mobility. On the other hand, results with 
Matern-cluster process showed that the estimation accuracy is poorer with clustered SBS 
deployment when compared to the homogeneous SBS deployment. Future research direc¬ 
tions may include: path prediction of UEs, and development of algorithms to dynamically 
adjust the measurement time interval depending on the past samples of estimated velocity. 


Appendix A Approximating the a and p Parameters 

of Gamma Distribution 

In this appendix, we derive the expressions for a and (3 parameters through a heuristic 
approach to minimize the MSE between the approximation f%{h) and the PMF fnih) of 
the handover count. First, we briefly describe a property of PPP with respect to the scaling 
of simulation space which will be used for the derivation of a and f3 parameters. 

Consider the square shaped simulation space shown in Fig. [T7|a) with an area L^. The 
SBS locations are shown in red dots, and the UE trajectory is shown as blue line of length 
d = vT. The intensity of SBSs can be expressed as A = where A^kvg is the mean 

number of SBSs in the simulation area. Assume that the simulation space is scaled down 
by the UE travel distance d, i.e., all the coordinates are scaled down by d, for example, 
coordinates of a SBS location (a;*, Ui) will be scaled down to (^, as shown in Fig. [T^b). 
Therefore, the simulation area will also be scaled down to (L/d)^. However, the mean 
number of SBSs A^avg remains the same, and hence the intensity of SBSs can be expressed as 
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(b) 

Figure 17: Illustration for scaling down the simulation space. 

A' = /L^ = XcP. Here, the term A' is a function of A and d, therefore, the statistics of 

the handover count H can be expressed in terms of A' alone. For example, the mean number 
of handovers in Q can be expressed as E[H] = dA/V/vr. Similarly, the parameters a and (3 
can also be approximated in terms of A'. 

We obtained histogram of the handover count H through extensive simulations in Matlab 
with a range of values for the SBS density A = {100,200,300,..., 10000} SBSs/km^ and the 
distance d = {0.033,0.1,0.2,0.4} km. For each combination of A and d, we obtained one 
million samples of H for constructing the PMF fnih). Then, we used the curve htting tool 
in Matlab to obtain the a and (3 parameters that provides a best £t of f%{h) in (|^ to the 
PMF fnih). By studying the characteristics of a and (3 parameters with respect to the 
variations in A and d, we formulated 


a = 2.7 + 4^)7 = 2.7 + dd^A, 

0.8 0.8 


(3 = 71 + 


= 71 + 


(34) 

(35) 


0.38+ \/V 0.38 + <i\/A 

The accuracy of these approximations are justified in Fig. nowhere the theoretical expressions 


(34) and (35) are shown to tightly overlap with the plots obtained through simulations. 
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Figure 18: Approximations of the a and (3 parameters. 


Appendix B Approximating the /i and Parameters 

of Gaussian Distribution 

In this appendix, we derive the expressions for /i and in ([^ to minimize the MSE between 
the approximation fnih) and the PMF fnih) of the handover count. We equate /r to the 
mean of handover count expressed in Q as 

AvTVX 4:d^/\ 

fj, = -=-. (36) 

TT TT 


Next, we obtain an expression for using a method similar to the one explained in Ap¬ 
pendix]^ We constructed the PMF fnih) through simulations, and then used curve htting 
tool in Matlab to obtain the parameter that provides a best fitting of in (|^ to the 

PMF fnih). By studying the characteristics of parameter with respect to the variations 
in A and d, we formulated 


a 


^ = 0.07 + 0.41 W = 0.07 + OAldVX. 


(37) 


The accuracy of this approximation is justified in Fig. 19 where the expression (37) is plotted 


in comparison with the plots obtained through simulations. 
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Figure 19: Approximation of the cx^ parameter. 


Appendix C Derivation of CRLB with Gamma Ap¬ 
proximation for Handover Count PMF 


In this appendix, CRLB for velocity estimation will be derived by considering the handover 
count PMF in (|^ that was derived using the gamma distribution. Taking the logarithm of 
PMF will provide the log-likelihood function, which will then be differentiated with respect 
to V. Finally, the derivative of log-likelihood function will be used to obtain the CRLB 
expression. Consider the handover count PMF, 


/&(*;«) 


r(a, /3/i, /?(/^ + 1)) 

r(5) 


(38) 


where, the a and j3 parameters are given by 0 and 0, respectively. By taking the logarithm 


of (38) and differentiating with respect to u, we get 


A 

dv 


log/i(h;u) 


A 

dv 


\ogT(a, I3h, I3{h + 1)) 


A 

dv 


logF(a). 


(39) 


Consider the hrst term in the RHS of (39). Let, zi 


jSh and Z 2 = (3{h + 1). Then, we 
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have 


A 

dv 


logr(a,/3h,/3(h + 1)) = 


d 


r(a, Zi, Z 2 ) dv 


r(a,Zi,Z2), 


d da d dzi 

N 7^r(a,Zi,Z2)^ + 7^r(a,zi,2;2) — 
i{a^zi^Z 2 ) \_da dw oz\ dv 


+ 


d , .dzs 

^r(a. 


(40) 


Each of the differentials in (40) can be derived to be 


A 

da 


T{a,zi,Z2) 


2 -f 2 (a, a; a + 1, a + 1; — Zi)zi 2 F 2 {a, a; a + 1, a + 1; — Z 2 )z 2 
a^ a^ 


- 7 (a, Zi) log( 2 :i) + 7 ( 0 , Z 2 ) log(z 2 ), 
^r(a, zi, Z2) = -e~^^z'j‘~\ ^r(a, zi, Z2) = e~^^Z 2 ~\ 

UZ\ OZ 2 

da — O.ShTo/X d ^2 —0.8(h + 1)T-\/A 

dv “ ’ dv “ (0.38 + vT^)2’ dv “ (0.38 + vT\/A)2 ’ 


where, 7 (a, x) = *df is the lower incomplete gamma function, and 2 -F 2 (ai, ^ 2 ; &i, ^ 2 ; z) 

is generalized hypergeometric function which can be expressed as 


2^2(01,02; fei, & 2 ; z ) — 

k=0 


{ai)k{a2)k z^ 
{bi)k{b2)k k\ ’ 


(41) 


where, (a)o = 1 and (a)^ = a{a + l)(a + 2)...(a + k — 1), for k > 1. 

Now, consider the second term in RHS of (39), ^logr(a) = It is well 

known that the logarithmic derivative of the gamma function is the digamma function '0(')- 
Therefore, 


g^logr(a)=V-(a)^ 




dv 
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Using the equations (40)-(42), we can express (39) as: 


^\ogffj{h-v) 


4TyA/3" 


a^V{a, I3h, I3{h + 1 )) 

4:T^A 

r (a, I3h, (3{h + 1)) 


a; a + 1, a + 1; —(3h) 

— (h + l)'^2-f2+ l,tt + 1; —/3(h 4“ 1)) 

-f{a,/3h) \og{/3h) - j{a,/3{h + 1)) log {/3{h + 1)) 

(43) 


0.8T^/XP'^-^e-^^ \h^ - e-^ih + 1)“1 ^ ^ ^ 

—--L^ — 4Ta/A 'ip{a). 

r(a,/4h,/4(h +l))(O.38 + nT0()2 


By squaring (43) and evaluating the expectation over H, we obtain the Fisher information 


as 


I{v) = E 


f d\og ffj{H;v) 
V dv 


(44) 


By taking the reciprocal of Fisher information, we can express the CRLB for h as in (12). 
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